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I. 


INTRODUCTION 


Several  independent  tasks  were  performed  under  this 
contract,  all  of  them  related  to  the  BRL-COMGEOM  technique o 
This  computer  technique  for  geometry  representation  was  orig- 
inally embodied  in  the  BRL-MAGIC  code  (written  by  MAGI)  and 
is  now,  in  a form  developed  by  BRL,  the  basis  of  the  BRL-GIFT 
code. 


The  tasks  completed  under  this  contract  and  described  in 
the  report  are: 

Ao  The  addition  to  the  list  of  available  body 
types  of  two  new  bodies,  the  general  ellip- 
soid and  the  "arbitrary  quadric  surface"? 

B.  The  development  of  a geometry  preprocessor 
program  capable  of  producing  GIFT  input 
from  a more  limited,  user-oriented,  set  of 
descriptions; 

C.  Development  of  a BRL-COMGEOM  model  of  the 
M109A1  tank. 

The  first  two  tasks  will  be  described  in  some  details 
The  modeling  task  is  embodied  in  an  input  deck  for  the  GIFT 
code  and  will  be  described  briefly  in  the  last  section. 
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II « ADDITIONAL  BODY  TYPES 


The  body  list  of  the  BRL-COMGEOM  (combinatorial  geoinetry) 
package  was  extended  to  include  the  general  ellipsoid  and  a 
quadratic  surface  generated  from  input  points » The  coding  that 
was  produced  for  each  of  these  bodies  covers  input  processing, 
ray  tracing  and  generation  of  unit  normals. 

Ao  The  General  Ellipsoid  (ELL) 

This  body  type  is  a generalization  of  the  ellipsoid  of 
revolution  which  was  included  in  the  body  list  of  the  original 
BRL-COMGEOM  package.  In  the  following  paragraphs  the  equations 
for  tracking  and  for  generating  the  local  surface  normal  will 
be  developed  and  the  code  implementation  will  be  described.  De 
tails  of  the  input  requirements  will  also  be  specified. 

1.  A Vector  Representation  of  the  General 
Ellipsoid. 

Let  X be  a point  on  the  surface  of  a general 
ellipsoid,  let  the  vertex  be  called  v,  and  let  the 
three  semi-axes  be  given  by 


a 

b 

c 


= aa^ 

= bb. 


= cc. 


(1) 


as  in  Figure  3. 

Now  we  define  three 

ing  x: 

a = (X  - v)  *aQ/a 
3 = (X  - v)*bQ/b 
Y = (x  - v)*Cq/c 


parameters  for  represent- 
(2) 
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FIGURE  1 


FIGURE  2 


FIGURE  3 


Ellipsoid  of  revolution, _specified  by  blank 
in  column  7 and  by  data  F^^,  ^2'^" 


Ellipsoid  of  revolution, _specified  by  1 in 
column  7 and  by  data  V,  a,  |5|, 


General  elli£soid,_specif ied  by  G in  column  7 
and  by  data  V,  a,  b,  c. 
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These  definitions  head  to  the  relationships; 

2.-2.  2 _ , 
a + B + Y 1 

ada  + 6d3  + Y<^Y  +0  (3) 

Thus  X can  be  expressed  in  terms  of  V,  a,  3,  Y 
and  a,  £ and  c: 

X = V + ,oa  + 3b  + YC  (4) 

2o  Ray  Tracing  to  the  General  Ellipsoid 

Let  a ray  emanate  from  a position  53  in  the 
direction  WBo_  We  wish  to  determine  the  point (s)  of 
intersection  X of  the  ray  with  the  general  ellipsoid 
and  the  distance  (s)  S from  XB  to  X,  We  define 

5 = 53  - V (5) 

and  note  that 

X - V = D + WB*S  . (6) 

Using  equations  (2)  for  a, 3 and  y we  find; 

{(D  + WB*S) =aQ/a}^  + {(D  + Vra»S)sbQ/b}^ 

+ {(D  + WB«S)»Cq/c)^  = 1, 

(7) 

Equation  (7)  is  a quadratic  in  the  distance 
S.  The  following  definitions  are  helpful  in  writing 
the  coefficients  of  the  quadratic; 

Cj^  = D-a^/a  ; C2  = V3*aQ/a 

c^  = D*bQ/b  ; c^  = WB»bQ/b 

Cg  = D^Cq/c  ; Cg  = WB»Cp/c  (8) 

In  terms  of  (8)  the  coefficients  of  the 
quadratic  term,  the  coefficient  of  the  linear  term 
and  the  constant  are,  respectively; 

2 2 2 

AA  = C2  + + Cg 

BB  = 2{C^C2  + + CgCg); 

CC  = - loO  o (9) 
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If  the  discriminant  is  less  than  zero  or 
vanishes,  the  ray  misses  the  ellipsoid  entirely  or 
just  grazes  it.  Two  meaningful  solutions  are  obtained 
if  the  discriminant  is  positive.  The  smaller  value 
is  identified  with  RIN,  the  distance  to  enter,  and 
the  larger  with  ROUT,  the  distance  to  emerge  from 
the  ellipsoid. 

The  numerical  technique  implemented  in 
the  ray  tracing  is  one  which  avoids  small  differences 
of  large  numbers  in  the  denominator. 

The  general  solution 


S 


(- 


can  also  be  written 


4AA‘CC  )/2AA 


2CC/ (-BB 


4AA*CC  ) . 


(10) 


(11) 


The  choice  between  (10)  and  (11)  is  made 
by  examining  the  sign  of  BB.  Let 


U 

U 


■J' 


BB  + /BB  - 4AA»CC 
- BB 


-J- 


BB  “ 4AA*CC 


BB<0 

BB>0 


Then 


S^=  U/(2AA) 
82=  2CC/U 


(12) 


is  a numerically  sound  solution  pair. 

3.  The  Unit  Normal  for  the  General  Ellipsoid 

The  unit  normal  to  the  ellipsoid  surface 
at  the  point  X can  be  obtained  as  the  normalized  cross 
product  of  two  differential  vectors,  both  tangent  to 
the  surface  at  X.  The  general  tangent  vector-  can  be 
obtained  from  (4)  and  from  (3) ; 
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(13) 


= ada  + bdg  - c (ada/y  + 3d3/y)  • 

If  we  choose  dX,  such  that  d3=0  and  dX, 
such  that  dtt“0,  then  ^ 

dXj^  = da  (a  - ac/y) 

_ _ _ (14) 

dX2  = d3(b  - 3c/y)  , 

and  the  unit  normal,  without  regard  to  sign,  is 
given  by 

n = dX^  X dX2/  dX^^  X dX2 

= (a’i^/a  + 3*bQ/b  + yc^/c) 


+ 3^b^  + yVc^  o (15) 

4.  Implementation  of  the  General  Ellipsoid 

The  coding  provided  for  the  implementation  of 
the  general  ellipsoid  consisted  of  an  input  routine  ELLP 
a tracking  routine  ELL  and  a section  of  coding  in  the 
routine  for  computing  normals. 

Three  properly  distinguished  forms  of  input 
are  acceptable  to  the  subroutine  ELLP.  Two  of  these 
describe  the  ellipsoid  of  revolution:  1)  the  two-focus, 

length-of-string  form  distinguished  on  the  leading  data 
card  by  a blank  in  column  7 and  2)  the  vertex,  semi- 
major axis  (vector) , semi-minor  axis  (magnitude)  form, 
distinguished  on  the  lead  data  card  by  any  entry  other 
than  a blank  or  a "G"  in  column  7o  The  third  form  is 
the  general  form  and  is  given  by  the  vertex  and  the 
three  (vector)  semi-axes.  The  general  form  is  signaled 
on  the  lead  data  card  by  a "G"  in  column  7.  The 
various  input  forms  are  demonstrated  in  Figures  1-3, 

The  function  of  ELLP  is  to  recognize  the  data 
form  and,  if  the  form  is  not  general,  to  convert  to  the 
general  form.  The  data  stored  by  ELLP  in  the  ASTER 
array  are  the  vertex  and  the  three  quantities  a^/a, 

bp/b,  Cq/c,  ELLP  also  stores  a pointer  for  locating 

these  data. 


Ray  tracing  as  performed  by  the  subroutine 
ELL  consists  in  straightforward  implementation  of 
equations  (5)  through  (12).  A miss  is  rejcorded  by 
returning  the  pair  RIN=+PINF  and  ROUT=-PINF,  where 
PINF  is  a very  large  number.  For  a hit  the  solution 
pair  is  assigned  so  that  ROUT  is  greater  than  RIN. 

B,  The  Arbitrary  Quadric  Surface 

This  body,  or  surface,  is  a completely  new  type  of  geo- 
metrical unit  and  requires  special  input  and  tracking  techniques 
which  are  discussed  below.  It  is  intended  to  permit  modeling 
of  elements  -which  are  not  approximated  by  any  of  the  other  body 
types,  but  which  can  be  fitted  in  sections  with  paraboloids. 

The  input  for  each  section  is  a set  of  measured  points  to  which 
a least  squares  fit  is  to  be  applied.. 

The  choice  of  a paraboloid,  as  opposed  to  a general  quadra- 
tic, is  related  to  the  region  description  technique:  we  need 

only  a small  section  of  the  surface  produced  but  we  must  insure 
that  the  rest  of  the  surface  does  not  re-enter  the  region  being 
described. 


1,  General  Considerations 

As  indicated  above,  the  technique  adopted 
here  is  not  a completely  general  quadratic  fit  but 
one  which  anticipates  the  function  to  be  performed 
by  the  surface  description  and  also  the  uncertainty 
attendant  on  any  fitting  technique.  This  is  illus- 
trated in  Figure  4,  in  which  the  AQS  forms  one 
boundary  of  a region  bounded  elsewhere  by  an  ellip- 
soid, The  quadratic  of  Figure  4a  doesn't  re-enter 
the  ellipsoid  and  the  region  description  is  there- 
fore meaningful.  However,  a general  quadratic 
could  result  in  a flat  ellipsoid  as  in  Figure  4b. 

Then  the  region  description  R=F+G  no  longer  des- 
cribes the  region  of  interest. 

The  paraboloid,  which  doesn't  close  on 
itself,  represents  our  best  means  of  avoiding  an 
ambiguous  region  description.  We  can  further  in- 
sure against  this  catastrophe  by  creating  the  para- 
boloid 2=h(u,v)  in  a co-ordinate  system  u,v,h  such 
that  the  u-v  plane  is  the  best  plane  through  the 
set  of  defining  points  - in  the  least  square  sense. 

This  approach  does,  however,  lead  to  a re- 
striction on  the  use  of  the  AQS.  A section  of  sur- 
face to  be  fitted  should  meet  the  requirement  that 
the  maximum  amount  by  which  the  section  deviates  from 
a plane  is  small  compared  with  the  radius  of  curva- 
ture of  the  section.  As  an  example  of  this  restric- 
tion, a paraboloidal  fit  to  a section  of  a sphere 
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REGION  R = 


FIGURE  4a:  AQS  (BODY  G)  does  not  reenter  BODY  F. 


REGION  R 


FIGURE  4b:  AQS  (BODY  G)  reenters  BODY  F and  truncates 

REGION  R. 


defined  by  a 45  degree  polar  angle  would  deviate 
from  the  sphere  by  approximately  0.05  of  the  radius 
at  the  apex  if  the  edge  is  matched  exactly.  , The 
ratio  of  the  maximum  height  of  this  section  to  the 
radius  (of  curvature)  is  0.3,  so  the  error  in  fitting 
is  about  16%  at  the  apex.  In  practice,  the  error 
would  be  divided  between  the  apex  and  the  edges. 

If  the  ratio  mentioned  above  is  limited  to  0.13 
(corresponding  to  a cone  of  30®  half  angle)  the 
error  in  fitting  is  about  a half  of  one  percent. 

2.  Theory  and  Analysis 

The  general  approach  is  to  assume  that  the 
piece  of  surface  represented  by  the  data 
{(X  ,Y  , Z ),  P=l,2,oooNP}  is  not  grossly  different 

r r 

from  a plane,  and  to  seek  first  the  "best"  plane, 
i.e.  that  plane  from  which  the  mean  square  perpen- 
dicular distance  of  the  points  is  a minimum.  This  is 
fruitful  in  several  ways:  if,  for  each  point,  the 

perpendicular  distance  from  the  plane  is  less  than 
the  input  tolerance,  T,  the  plane  can  be  accepted  as 
a proper  representation.  If,  however,  the  plane  is  not 
sufficiently  representative,  it  can  be  used  as  one 
of  the  co-ordinate  planes  of  a co-ordinate  system 
in  which  the  paraboloidal  representation  of  the 
surface  is  single-valued: 

2 2 

h(u,v)  = aj^u  + a2^  + a^u 

+ a^v  + Ug  (16) 

where  the  u-v  plane  is  the  "best"  plane  and  the 
third  co-ordinate  h is  the  actual  height  of  the 
surface  above  the  best  plane.  The  difference 
{h (u^,v^) -h^}  for  the  individual  point  will  not  be 

very  different  from  the  perpendicular  distance  of 
the  point  from  the  surface  and  can  be  used  as  a 
measure  of  the  success  of  the  fit. 

The  first  step  in  processing  is  to  perform 
a transformation  into  the  center  of  mass  co-ordinate 
system.  This  is  simply  a computational  convenience 
and  involves  subtracting  from  each  point  vector, 
rp,  the  mean  value  of  rp  over  all  input  points.  In 
this  system  the  best  plane  is  sought. 
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The  method  is  to  minimize  the  mean  square 
perpendicular  distance  of  points  to  the  plane 

subject  to  the  condition  that  3^^  + ^2^  + 3^^  = 1 (18) 

or  that  these  3's  are  direction  cosines. 

The  perpendicular  distance  of  a point 
from  the  plane  is 

. + (19) 

and 

ZEr  = 2 (Bj^Xr  + B2Yk 

Since  there  is  no  constraint  involving  34  we  can 
solve  directly  for  it: 

2 

9Zek 

= 0 = 2(Bj_  ZXr  + B2  SyR  + B3  ZZr  + 34'NP)  (21) 

In  the  center  of  mass  system  the  sums  on  the 
right  hand  side  of  equation  (21)  all  vanish  - whence 
= 0 and  the  plane  passes  through  the  center  of 

mass  (as  one  might  have  guessed  to  start  with) . • 

Let  us  call  the  i co-ordinate  of  the  k^ 
point  w^R  and  write  that 


"k  = I ei'^ik 


1 T-  2 

!?P 

k 


- e - 


NP  J ? ^^i^j'^ik'^jk 


Z ZBiA,jBj 


(22) 


(23) 
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(24) 


Where  the  E w.,^w.)^  = w^w^ 


are  the  elements  of  a symmetric  matrix  A,  Then  3 can 
be  thought  of  as  a vector  in  3-space  and 


= (BA)B. 


(25) 


The  minimization  of  e'^  subject  to  the  con- 
2 

strarnt  B =1  requires  that 


3e^  _ ,3B^  _ 0 


for  £ = 1,2,3. 


(26) 


or 


2Z  A,.Bj 


- 2ABj,  = 0 


(27) 


or 


(A-AI)B  = 0 


(28) 


where  I is  the  unit  matrix. 


A is  then  the  solution  of  a cubic  equation. 
Its  physical  significance  is  seen  by  multiplying 
the  condition  equation  by  B^^  and  summing  over  £ : 


or 


AZB.^  = 0 = (BA)B-A 
£ ^ 


(29) 


Thus  the  three  values  of  A correspond  to 

2 

the  local  extrema  of  the  quantity  e under  rotation 
of  the  plane.  The  smallest  value  of  A gives  the  de- 
sired plane. 

After  the  plane  has  been  determined  the 
action  of  the  processing  subroutine  AQSIN  is  to  com- 
pare the  perpendicular  distance  of  each  point  with 
the  tolerance  value  to  determine  whether  the  plane 
should  be  accepted.  If  it  is  accepted,  the  quadric 
fit  section  described  in  the  next  few  paragraphs  is 
by-passed  and  the  coefficients  are  corrected  to 
refer  to  the  laboratory  coordinates. 
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If  the  plane  is  not  accepted  the  para- 
boloidal fit  is  pursued.  A transformation  is  made 
(a  rotation)  to  a coordinate  system  u,v,h  such  that 
the  "best"  plane  is  the  u-v  plane  and  h{u,v)  will 
describe  the  paraboloid: 

2 2 

h{u,v)  = ctj^u  +0.2^  + 0^4^ 

+ aeV  + a6  = Za  .w.  ; = 1 (30) 

5 j ^ ^ ° 

The  quantity  to  be  minimized  is  the  mean 
square  of  the  point  height  above  (or  below)  the 
paraboloid.  Under  the  assumed  conditions  this  differs 
little  from  a perpendicular  distance  and  is  much 
easier  to  handle: 


\ (31) 

- h2]  <32) 


Again  we  identify  the  symmetric  matrix  A whose 
ij  element,  is 


^ij  N ^ '^ik'^jk  ; 
and  the  vectors  a and  w such  that 


= (aA)a  - 2h(aw)  + h^  ; 


and 


= O' 


(33) 


(34) 


(35) 


In  vector  notation: 
a = A ^ (hw) 


(36) 
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Thus  the  code  constructs  the  matrix  A: 


(37) 


(38) 


The  matrix  equation  is  solved  by  the  sub- 
routine SOLVE  to  produce  the  6-element  vector  a. 


The  vector  a consists  of  the  coefficients 
of  the  paraboloid  in  the  rotated  center  of  mass 
system.  The  inverse  transformations  to  the  labora- 
tory system  must  be  performed  to  produce  the  10 
element  vector  of  coefficients  in  that  system. 


The  10  coefficients  of  the  equation 


Bq  + Bj^x  + B2y  + B^z  + B^xy  + B^xz  + B^yz 

+ B^x^  + Bgy^  + BgZ^  = 0 (39) 

are  the  essential  data  stored  for  the  representa- 
tion of  the  AQSo 


13 


3.  Ray  Tracing  to  the  AQS 

Ray  tracing  to  the  AQS  consists  in  ^ substi- 
tuting into  (39)  the  following  expressions  for  the 
coordinates  of  the  point  of  intersection; 

X = + WB«S  , (40) 

where  (40)  is  analogous  to  (6)  for  the  general 
ellipsoid.  This  substitution  results  in  a quadratic 
in  S which  is  solved  in  a manner  entirely  analogous 
to  that  described  for  the  general  ellipsoid  (equations 
10  through  12). 

The  roots  of  this  quadratic,  and  82^  must 

be  interpreted  with  care  because  the  AQS  is  not  a 
closed  convex  body.  Figure  5 demonstrates  the  various 
special  cases  to  be  considered  and  the  interpretation 
of  and  Sj  (S2>S^)  as  RIN  cind  ROUT. 

a)  An  AQS  which  is  concave  toward  the  inside 
has  a single  pair  specified  for  the  inter- 
sections; (RIN, ROUT)  = (Sj^,  S2) 

b)  An  AQS  which  is  convex  toward  the  inside 
has  a pair  associated  with  each  intersection: 

(RIN, ROUT)  = (-PINF,Sj^)  and 
(RIN, ROUT)  = (S2/PINF) 

c)  A plane  AQS  is  covered  by  (a) ; 

Entering;  (RIN,  ROUT)  = (Sj^,PINF) 

Leaving;  (RIN, ROUT)  = (-PINF,S2) 

4.  The  Normal  to  the  AQS 

The  normal  calculation  is  straight  forward 
and  consists  in  taking  the  gradient  of  the  function 
given  by  equation  (39); 

Vf(x,y,z)  = 1 (Bj^  + B^y  + B^z  + B^x) 

+ 5 (B2  + B^x  + BgZ  + Bgy) 

+ ic(B3  + BgX  + Bgy  + BgZ)  (41) 

and  normalizing  it.  The  signs  of  the  coefficients  in 
(39)  have  been  chosen  so  that  the  normal  points  toward 
the  "inside",  as  indicated  by.  the  point  given  on  input. 
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ray 


a.  A concave  AQS,  (RIN,ROUT)  = (Sj^/S2) 


b,  A convex  AQS, 

(RIN,ROUT)j^  = (-PINF,S^) 
(RIN,R0UT)2  = (S2,PINF) 

FIGURE  5:  Tracking  to  the  AQS 
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5.  Implementation 

Three  subroutines  were  developed  to  handle 
the  AQS  in  the  GIFT  code  and  two  old  routines  were 
altered: 

a.  AQSIN(LBOT,  DATAjFX)  is  a new 
subroutine  called  by  GENI  to  read 
and  process  AQS  data. 

b.  SOLVE (A,B,KS,N)  is  an  IBM  sub- 
routine called  by  AQSIN  to  invert 
a matrix  and  return  coefficients 
of  the  quadratic  in  a special  co- 
ordinate system, 

c.  AQS  was  added  to  implement  the 
solution  of  the  ray  tracing  equation. 

d.  GENI  was  altered  to  accept  the  AQS 
input  and  to  call  AQSIN. 

e.  CALC  was  altered  to  compute  the 
normal  for  the  AQS, 

The  input  for  the  AQS  consists  of: 


Card  1/ 

read  by  GENI  in  format  (3A1,  A3, 

A4,  6F10.0) 

Columns 

1-3 

(body  number) 

4-6 

AQS 

11-20 

NP,  number  of  points  (fixed  in 

AQSIN) 

21-30 

TOL,  tolerance 

31-40 

41-50 

Y>  coordinates  of  an  "inside 

" point 

51-60 

Z ) 

Subsequent  cards  are  read  by  AQSIN  in  format 

(lOX,  6F10.5)  and  contain  coordinates  of 

pairs  of  points 

Columns 

11-20 

Xp 

21-30 

Y 

P 

31-40 

Z 

p 

41-50 

Sr 

Vi 

51-60 

Vi 

61-70 

^p+1 
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The  points  are  read  into  the  ASTER  array  in 
which  6*{NP+8)  temporary  locations  are  reserved  for 
processing. 

The  data  stored  in  the  ASTER  array  for  the 
AQS  (ITYPE=12)  of  body  number  NBO  is  the  following: 


LOCATION 

LBODy+3* (NRPP+NBO-1) 
LDATA 

LDATA+1 

LDATA+2 

LOCDA 


CONTENTS 

2^^*ITYPE+LDATA 

Location  of  first  coefficient 
LOCDA 

2^^* (LOCDA+10)+ (LOCDA+13) 

2^^* (L0CDA+16)+ (LOCDA+19) 

. 2 
coefficient  of  x 


LOCDA+9 


Bq/  constant  term 
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III.  BRITL:  A GEOMETRY  PREPROCESSOR  PROGRAM 


FOR  INPUT  TO  THE  GIFT  SYSTEM 

A.  Introduction 

BRITL  (Body-Region  Input  Translator)  was  developed  as  an 
aid  in  preparing  geometry  data  for  the  GIFT  program.  It's 
several  functions  are  by  no  means  exhaustive  of  the  possibili- 
ties but  represent  an  important  beginning  in  the  automation  of 
input  production. 

The  tasks  of  the  preprocessor  are  listed  below  and  are 
discussed  in  some  detail  in  the  subsequent  sections: 

1.  Given  a body  (other  than  an  ARS)  from  among 
the  types  acceptable  to  the  GIFT  program,  BRITL  produces  an 
inner  body  of  the  same  type  by  moving  the  surfaces  (nearly) 
normal  to  themselves  by  distances  given  by  input.  Both  bodies 
are  recorded  in  the  body  table  for  use  together  as  a "thick" 
body. 


2.  Given  a series  of  vector  positions  and  associated 
radii,  BRITL  forms  a PIPE  by  joining  successive  positions  with 
either  an  RCC  or  a TRC,  as  indicated  by  the  comparison  of  radii. 

A PIPE  may  be  thickened. 

3.  BRITL  produces  region  descriptions  (as  called  for 
in  the  region  table)  corresponding  to  the  "shell"  of  the  "thick" 
body,  the  interior  or  the  exterior  of  the  "thick"  body. 

4.  BRITL  allows  backward  references  to  already  de- 
fined regions  and,  by  applying  DeMorgan's  rules,  produces  a 
region  description  in  the  required  conjunctive  form. 

5.  Given  a series  of  regions,  already  described, 

BRITL  produces  a new  series  of  regions  by  means  of  a transforma- 
tion applied  to  each  of  the  bodies  in  the  description.  In 
general,  the  transformation  must  consist  of  rotations,  trans- 
lations and  reflections  without  distortion  or  scale  change. 

For  special  bodies  (those  described  entirely  by  vectors)  dis- 
tortion is  allowed. 

6.  BRITL  allows  the  stripping  of  "model"  regions  and 
unused  bodies  from  the  body  and  region  tables.  It  produces 
properly  sequenced  and  numbered  bodies  and  regions  and  region 
identifications  in  the  format  required  by  the  GIFT  program. 

Body  and  region  numbers  may  start  at  initial  values  provided  on 
input. 

Because  new  bodies  are  produced,  and  because  region  des- 
criptions may  refer  to  previous  regions  as  well  as  to  resequenced 
bodies,  numerical  sequencing  of  bodies  and  regions  is  not  feasible. 
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Instead  bodies  and  regions  are  "neimes"  with  distinct  3-character 
alphanumeric  tags.  The  reference,  within  a region  description, 
to  a body  or  region  is  made  with  a signed  four-character  name,  the 
last  three  characters  of  which  is  the  tag  mentioned  above » The 
first  character  is  used  in  connection  with  the  "thick"  body  to 
designate  shell,  interior  or  exterior. 

Thus,  the  input  to  the  geometry  preprocessor  is  1)  a body 
table,  with  bodies  designated  by  names,  and  with  some  bodies 
designated  as  "thick";  2)  a region  table  with  references  to  bodies 
and  previously  described  regions,  and  special  regions  designated 
as  transformations  of  previously  described  regions  and  accompanied 
by  transformation  data  and  3)  a region  identification  table, 
which,  in  addition  to  the  usual  information,  may  carry  a "model" 
tag.  The  input  formats  are  as  close  as  possible  to  those  of  the 
GIFT  program  and  the  output  formats  are  precisely  those  of  the 
GIFT  program. 

The  basic  formats  are  given  here: 

Title  Card  (20A4) 

Parameter  Card  (3110) 

Col\nnn  1-10  Number  of  enclosing  RPP*s 

11-20  Number  of  bodies 

21-30  Nvimber  of  regions 

Enclosing  RPP  card  (2A3,  A4,  6E10.3,  A4,  2A3) 

Coliimns  1-10  Comments 

71-80  Comments 

11-70  Standard  RPP  data 

Lead  Body  Card  (2A3,  4A1,  6F10.0,  A4,  2A3) 

Columns  1-3  Body  name  (alphanumeric) 

Columns  4-6  Body  type  (ARB,  REC,  etc.) 

7 Used  to  distinguish  ARB  types, 

otherwise  blank 
8-9  Blanks 

10  T for  a thick  body,  otherwise  blank 

11-70  Standard  body  data 

71-80  Comments 

Succeeding  Body  Cards  (lOX,  6F10.5,  A4,  2A3) 

Same  as  lead  body  card  except  that  first  ten  columns 
are  not  read. 

Thickness  Card(s)  (lOX,  6F10.5,  A4,  2A3) 

Same  format  as  body  cards  - six  thicknesses  per  card. 

Read  only  if  T in  column  10  of  lead  body  card. 

Region  Card  (2X,  A3,  IX,  9(A2,  2A1,  A3)) 

Columns  3-5  Region  name  (alphanumeric) 
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For  the  I-  body  or  region  reference,  first  column 
is  5+7.  (I-D+l  = NI  + 1 

Columns  (NI+1)- (NI+2)  OR  or  blanks 

Column  (NI+3)  Sign  (+  or  -) 

Column  (NI+4)  Special  symbol  for 

reference  to  a "thick"  body 

Columns  (NI+5) - (NI+7)  3 character  body  or  region 

name 

For  the  special  case  of  a region  or  series  of  regions,  to 
be  obtained  by  transformation  from  a previous  region 
series,  the  region  card  has  the  above  format  for  the  follow- 
ing information: 

Columns  3-5  Region  name  (one  name  for 

whole  series) 

Columns  7-8  DP 

Columns  11-13  Name  of  first  region  in 

the  model  set 

For  this  case  the  region  card  is  followed  by  a card 
bearing  the  number  of  regions  (15)  and  by  a set  of  four 
cards  with  transformation  data  for  four  non-coplanar 
points,  (lOX,  6F10.5): 

Columns  11-20  X ^ 

21-30  Y ^ point  in  space  of  model  regions 

31-40  Z 

41-50  X' 

51-60  Y' 

61-70  Z' 

End  Card  for  Region  Table;  -1  in  colvunns  4-5 

Region  Identification  Cards 

(A3,  A2,  2X,  A3,  A3,  A4,  A3,  IlO,  lOX,  10A4) 

Column  1-3  Either  blank 

or  MOD  for  a model  region 

or  DUP  for  a region  set 

obtained  by  trans- 
formation 

Columns  4-5  Either  blank 

or  ST  if  this  is  a region 
to  be  stripped  out 
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Columns  8-10 
Columns  11-20 


Region  name 

Ident  code  - right  justified 
or  for  a region 
set  obtained  by 
transformation , 
name  of  first 
region  that  is  to 
be  transformed  - 
right  justified 

Columns  21-30  Space  code  - right  justified 

or  for  a region 
set  obtained  by 
transformation , 
number  of  regions 
to  be  transformed 
right  justified 

Columns  41-80  Region  comments 

Initial  body  number  and  Initial  region  number  card  (2110) 

Columns  1-10  Integer  indicating  the  niomber 

of  the  first  output  body 

Columns  11-20  Integer  indicating  the  niomber 

of  the  first  output  region 

B.  The  Thick  Body 

Any  body  type  acceptable  to  the  GIFT  code,  except  the  ARS, 
can  be  "thickenedo " In  addition,  the  PIPE  body  can  be  "thickened. 
The  number  of  thicknesses  that  can  be  specified  for  a given  body 
type  is  peculiar  to  that  body  as  described  below,  A "T"  in 
column  7 of  the  lead  data  card  for  the  body  indicates  to  the  code 
that  the  body  is  to  be  thickened,  i.eo,  that  another  body  is  to 
be  produced.  Thickness  card(s)  follow  the  body  data  cards  (up 
to  6 thicknesses  per  card  in  the  same  format  as  the  body  cards) 
except  for  the  PIPE  body.  The  order  in  which  the  thicknesses 
are  given  associates  them  with  the  proper  dimensions o 

1.  BOX  (also  RPP) 

Six  distinct  thicknesses  can  be  specified,  one  for 
each  side,  in  the  order  shown  in  Figure  6a,  Thus  t^i  and  t^ 

have  the  effect  of  moving  the  vertex,  V,  along  each  of  the  height 
vectors  in  order  (and  of  shortening  the  height  vectors)  and  t2/ 
t^  and  tg  further  shorten  the  height  vectors  H2  and 

The  body  type  RPP  is  converted  to  a BOX  before  thickening. 
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2o  RAW 


The  right  angle  wedge  has  almost  the  same  thickness 
format  as  the  BOX  except  that  the  inputs  and  t^  are  equal  to 

each  other  and  to  the  distance  that  the  incl^ined  face  is  to  be 
moved  normal  to  itself <,  The  shortening  of  and  H2  are  calcu- 
lated internally.  See  Figure  6b. 

3.  SPH 

The  sphere  can  have  one  thickness  along  the  radius. 

4 o RCC 

The  RCC  can  have  three  thicknesses:  t^^  normal  to  the 
base,  t2  normal  to  the  upper  face,  t^  normal  to  the  curved 
surface o 

5.  REC 

The  REC  can  have  three  thicknesses,  just  similar  to 
those  of  the  RCC.  The  3rd  thickness  is  interpreted  as  a short- 
ening of  the  semi-axes. 

6.  TRC 

The  TRC  has  three  possible  thicknesses,  as  for  the 
RCC  and  REC. 

7.  ELL 

The  ELL  has  one  thickness,  tj^,  interpreted  as  a short- 
ening of  the  semi-axes. 

8 o ARB 

The  ARB  can  have  one  thickness  for  each  face,  given  in 
the  same  order  as  the  face  descriptions. 

There  is  one  pathological  case,  namely  that  of  the 
five-sided  pyramid.  The  problem  arises  because  the  apex  is  the 
intersection  of  four  planes  (rather  than  three) . A general  choice 
of  thicknesses  would  lead  to  one  of  several  different  ARB  types. 
This  complication  has  been  avoided.  It  has  been  ass\amed  that 
the  inner  body  is  congruent  with  the  original  and  the  intersection 
of  the  first  three  displaced  planes  is  taken  as  the  apex.  Thus 
the  only  meaningful  input  thicknesses  are  those  that  preserve 
the  shape . 

9.  TEC 

The  TEC  can  have  three  thicknesses,  described  as  for 
the  RCC,  REC  and  TRC. 
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10 . TOR 

The  TOR  has  one  thickness  corresponding  to  a short- 
ening of  the  small  radius. 

lie  ARS 

No  provision  has  been  made  for  thickening  this  body. 

12.  AQS 

The  AQS  may  have  either  one  thickness  (indicated  by 
a zero  in  columns  61-70  of  the  lead  card)  or  one  thickness  for 
each  input  point  (indicated  by  a 1.0  in  columns  61-70).  The 
action  of  the  code  is  to  form  an  AQS  through  the  original  points, 
calculate  a local  normal  (approximate)  for  each  point,  and  then 
create  a new  point  (belonging  to  the  new  AQS)  at  the  proper  dis- 
tance. Thus  the  output  AQS  is  just  a set  of  points  derived  from 
the  original  body  and  its  thickness  set. 

C.  Region  Descriptions  Based  on  the  Thick  Body 

Let  the  body  to  be  thickened  have  the  name  ABC.  A refer- 
ence to  the  shell  (see  Figure  7)  in  a region  description  is 
+SABCo  A reference  to  the  interior  is  +IABC. 

If  the  translator  gives  the  original  body  the  number  N the 
inner  body  will  have  the  number  N+1  and  SABC  will  be  replaced  by 
{N-(N+1)}  while  lABC  will  be  replaced  by  N+1  . 

The  former  references  are  handled  automatically  in  the  code, 
but  two  other  reg[ions  demonstrated  in  Figure  7b,  can  be  referred 
to  by  the  user  in  a region  description  in  various  ways; 

lo  For  the  region  BRG  between  ABC  and  an  enclosing 
body  BBC  (numbered  NBB  by  the  code),  the  best  (i.e.  least  re- 
dundant) input  description  is 

BRG  = BBC  - ABC 

which  translates  to  {NBB-N} 

2.  The  union  of  BRG  with  lABC  should  be  referenced  as 
OR (BCC-ABC ) OR ( lABC) 
which  translates  to  OR (NBB-N) OR (N+1) . 

The  other  possible  descriptions  of  these  regions  have 
geometrical  redundancies  which  are  not  recognized  by  the  trans- 
lator. 
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D,  The  Pipe 


The  pipe  consists  of  a series  of  truncated  cones  or  cylinders 
connected  by  spheres,  as  shown  in  Figure  8.  The  options  available 
for  the  pipe  are  as  follows: 

1,  constant  radius  (cylindrical  sections) 
or  variable  radius  (conical  sections); 

2,  walls  of  non-zero  thickness  (inner 
bodies  are  constructed  with  the  normal 
separation  constant  in  each  section)  or 
zero  thickness  (no  inner  bodies  needed) ; 

3,  for  a thick  pipe  the  wall  thickness  may 
be  the  scune  or  different  for  all  sections; 

4,  the  pipe  may  close  on  itself  or  it  may  be 
double-ended . 


The  essential  input  to  the  PIPE  is  a set  of  NP  vector  posi- 
tions, radii  and  (possibly)  thicknesses:  , p=l,2...NP) 

XT  XT  ST 

Additional  input  on  the  lead  data  card  gives  the  common  radius 
if  the  pipe  is  cylindrical,  and  the  common  thickness  if  the  pipe 
walls  have  the  same  normal  thickness  in  each  section.  Zero  values 
of  these  parameters  indicate  that  the  radius  and  thickness  are 
variable.  Details  of  the  input  are  given  at  the  end  of  this 
section. 


1)  SNAME  can  refer  only  to  a thick  pipe  and  means 
the  shell  of  the  pipe.  If  the  outer  bodies  are 
designated  as  and  the  inner  bodies  as 

Cj,  I=1,2...NBD  as  in  Figure  2,  the  region 

SNAME  is  the  union  of  the  regions: 

®I  “ " *^i+i 


for  all  I such  that  1 < I < NBD  and 


and 


C2  (1=1) 


®NBD  “ ^NBD-1  “ *^NBD  “ ^1 


(I=NBb) 


This  specification  handles  both  the  double- 
ended  pipe  and  the  closed  pipe. 


2)  INAME  is  a valid  reference  for  a thick  pipe  or 
a thin  pipe.  In  either  case  it  means  the  inner 
air  region  and  is  given  by  the  union  of  all  C ^ 
(or  of  all  Bj  for  a thin  pipe). 
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3)  NAME  is  a valid  reference  for  either  a thick 
pipe  or  a thin  pipe,  but  is  not  very  useful  for 

a thick  pipe.  It  is  identical  to  INAME  for  a thin 
pipe.  It  is  the  union  of  all 

4)  -NAME  is  valid  and  meaningful  for  either  a 
thick  pipe  or  a thin  pipe  and  refers  to  the  out- 
side region.  It  is  the  intersection  of  all 
negated 'Bj, 

The  action  of  the  preprocessor  program  is  the  following; 
the  lead  data  card  is  read  and  its  information  stored  in  the 
ASTER  array.  When  the  body  has  been  identified  as  a pipe,  the 
pipe  index  LPIPE  is  increased  by  1,  the  starting  body  number  is 
stored  in  the  array  NNPIPE (LPIPE)  and  the  subroutine  PIPE (LI, L2) 
is  called,  LI  is  the  starting  location  in  the  ASTER  array  of 
the  pipe  data  and  L2  will  be  returned  as  the  last  location. 

PIPE (LI, L2)  reads  the  remaining  data  into  temporary  storage,  con 
structs  the  bodies  (RCC's  or  TRC's  and  SPH's)  and  updates  all 
pointer  arrays  indexed  on  body  number.  In  addition,  it  counts 
the  new  bodies  and  stores  that  number  in  the  array  LBD (LPIPE). 

The  rest  of  the  coding  related  to  the  pipe  is  in  the  region 
processing  section  of  the  main  program.  When  a name  in  a region 
description  has  been  identified  as  a body  name,  the  body  number, 
NN,  is  compared  with  the  list  NNPIPE (LPIPE) . If  there  is  an 
LPIPE  such  that  NNPIPE  equals  NN,  control  is  transferred  to  the 
pipe  section  where  the  prefix  of  the  NAME  is  examined  to  deter- 
mine the  appropriate  region  description. 
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The  lead  data  card  for  a pipe  is  formatted  as  follows: 


COLUMN 

FORMAT 

1-3 

A3 

NAME  of  pipe  body  (alphanumeric) 

4-6 

A3 

"PIP",  identifier 

7-9 

3A1 

Blank 

10 

A1 

Blank  or  "T"  for  a thick  pipe 

11-20 

FIO.O 

Number  of  points 

21-30 

FIO.O 

Common  radius  for  all  points  or  0.0 

31-40 

FIO.O 

Common  wall  thickness  for  all  sections 
or  0.0 

41-50 

FIO.O 

0.0  for  double-ended  pipe,  any  number 
y 0.0  for  a closed  pipe 

51-70 

Blank 

71-80 

A4 , 2A3 

Comments 

The 

rest  of  the 

data  cards,  one  per  point,  specify  the  point 

data: 

COLUMN 

FORMAT 

1-10 

Not  read 

11-20 

F10.5 

21-30 

F10.5 

Y 

P 

31-40 

Flo. 5 

Z 

P 

41-50 

F10.5 

Radius-unless  a common  radius  has  been 
specified 

51-60 

F10.5 

Thickness  - unless  a common  thickness 
has  been  specified  or  the  body  is  thin 

61-70 

Blank 

71-80 

A4,2A3 

Comments 
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E,  Region  Transformations 


A series  of  regions  may  be  designated  as  "model"  regions 
and  a backward  reference  on  a region  card,  togetber  with  a 
card  giving  the  number  of  regions  in  the  series  and  four  cards 
giving  transformation  data  for  four  non-coplanar  points,  will  pro- 
duce a series  of  new  bodies  and  regions o The  format  for  these 
data  has  been  given  in  the  introduction.  Here  we  simply  show  a 
sample  input  (Figure  10)  and  the  resulting  output  (Figure  11) . 

A model  of  a man  (crew  member)  consists  of  11  bodies  and 
11  regions,  as  shown  in  Figure  10.  The  model  is  to  be  rotated 
and  translated  into  each  of  two  new  positions  (regions  RAP  and 
RAQ,  respectively) . Following  the  region  cards  for  RAP  and  RAQ 
is  a card  bearing  the  number  11  to  indicate  the  number  of  regions 
to  be  transformed.  The  four  following  cards  carry  the  4-point 
transformation  P " 1,2, 3, 4). 

The  region  identification  cards  for  the  model  bear  the  tag 
MOD  and  the  tag  ST  which  indicates  that  the  model  is  to  be 
stripped  out  at  the  end.  The  region  identification  cards  (one  per 
series)  for  the  copies  bear  the  tag  DUP,  the  name  of  the  first 
region  in  the  model  (RAA)  and  the  number  of  regions  in  the  model 
set.  Finally  there  is  a card  giving  the  initialization  of  body 
number  (536)  and  region  number  (536)  for  the  output. 

The  output.  Figure  11,  shows  properly  sequenced  body, 
region  and  region  ident  cards  for  the  two  crew  members  copied 
from  the  model.  The  model  itself  has  been  stripped  outo 

The  four-point  transformations  here  represent  translations 
and  rotations  only.  One  point  is  therefore  redundant  (and  must 
be  consistent).  For  regions  involving  only  the  AQS,  the  ARS, 
the  ELLG  and  the  TEC  scale  changes  and  distortion  are  allowed 
and  these,  in  general,  require  the  four  point  specification. 
Reflections  in  a plane  can  be  specified  for  any  body. 

The  output  of  BRITL  appears  on  TAPE  7 and  can  be  trans- 
ferred to  the  PUNCH  file  or  the  print  file  via  job  control  cards. 

F.  Subroutines  and  Important  Arrays 

In  this  section  we  list  and  describe  briefly  the  important 
arrays  and  the  subprograms  of  BRITL o 


1.  Some  important  arrays: 


ASTER-MASTER 

LOCDB(2000) 

LOCDR(IOOO) 

NAME (2000) 


currently  dimensioned  at  10000 

starting  locations  in  ASTER  of  body  data, 
indexed  on  body  number 

starting  location  in  ASTER  of  region  data, 
indexed  on  region  number 

array  of  body  names  (3  characters)  versus 
body  numbers 
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CREW  MEMBERS 

0 11  13 


AAAELL 

0.  0 

o 

• 

o 

16.5 

o 

. 

o 

o 

. 

o 

20.0 

HEAD  1 

9.0 

HEAD  2 

AABRCC 

0.0 

0.0 

9.4 

0.0 

0.0 

5.0 

HEAD  1 

2.0 

NECK  2 

AACREC 

0.0 

-5.75 

-0.425 

0.0 

11.5 

0.0 

BODY  1 

0.0 

0.0 

11.075 

4.0 

0.0 

0.0 

BODY  2 

AADRCC 

0.0 

-7.00 

-3.75 

0.0 

0.0 

14.4 

UPRARM  1 

1.25 

UPRARM  2 

AAERCC 

0.0 

7.00 

-3.75 

0.0 

0.0 

14.4 

UPLARM  1 

1.  25 

UPLARM  2 

AAFRCC 

1.25 

-7.00 

-2.50 

12.25 

0.0 

0.0 

LORARM  1 

1.25 

LORARM  2 

AAGRCC 

1.25 

7.00 

-2.50 

12.25 

0.0 

0.0 

LOLARM  1 

1.25 

LOLARM  2 

AAHRCC 

0.0 

-3.50 

-9.0 

16.5 

0.0 

0.0 

RTHIGH  1 

2.5 

RTHIGH  2 

AAIRCC 

0.0 

3.50 

-9.0 

16.5 

0.0 

0.0 

LTHIGH  1 

2.5 

LTHIGH  2 

AAJRCC 

15.0 

-3.50 

-9.0 

0.0 

0.0 

-15.0 

RT  LEG  1 

1.5 

RT  LEG  2 

AAKRCC 

15.0 

3.5 

-9.0 

0.0 

0.0 

-15.0 

LF  LEG  1 

1.5  LF  LEG  2 


FIGURE  10;  Sample  Input  - Body  Table 


AAA 


AAC 


OJ 

NJ 


RAA  AAA 

RAB  AAB  - 

RAC  AAC 

RAD  AAD  - AAC 

RAE  AAE  - AAC 

RAF  AAF  - AAC 

RAG  AAG  - AAC 

RAH  AAH  - AAC 

RAI  AAI  - AAC 

RAJ  AAJ  - AAH 

RAK  AAK  - AAI 

RAP  DP  RAA 

11 

0.0 

13.5 

13.5 

15.0 

RAQ  DP  RAA 
11 

0.0 

13.5 

13.5 

15.0 

-1 


MODST 

RAA 

10 

MODST 

RAB 

11 

MODST 

RAC 

12 

MODST 

RAD 

13 

MODST 

RAE 

14 

MODST 

RAF 

15 

MODST 

RAG 

16 

MODST 

RAH 

17 

MODST 

RAI 

18 

MODST 

RAJ 

19 

MODST 

RAK 

20 

DUP 

RAP 

RAA 

DUP 

RAQ 

RAA 

536 

536 

- AAD 

- AAE 


0.0 

-7.0 

7.0 

-3.5 


0.0 

-7.0 

7.0 

-3.5 


11 

11 


FIGURE  10,  continued 


HEAD 
NECK 
BODY 
UPRARM 
UPLARM 
LORARM 
LOLARM 
RT  THIGH 
LF  THIGH 
RT  LEG 
LF  LEG 


A 

A 

0.0 

15.206 

-52.0 

-12.85 

A 

-2.5 

22.206 

-38.5 

-15.35 

A 

-2.5 

8.206 

-38.5 

-15.35 

A 

-9.0 

18.708 

-37.0 

-21.85 

A 

B 

B 

0.0 

0.0 

-52.0 

-12.85 

B 

2.5 

7.0 

-38.5 

-15.35 

B 

2.5 

-7.0 

-38.5 

-15,35 

B 

-9.0 

3.5 

-37.0 

-21.85 

B 

HEAD 

NECK 

BODY 

UPPER  ARM  RIGHT 
UPPER  ARM  LEFT 
LOWER  ARM  RIGHT 
LOWER  ARM  LEFT 
THIGH  RIGHT 
THIGH  LEFT 
LOWER  LEG  RIGHT 
LOWER  LEG  LEFT 

CREW  A 
CREW  B 


Region  table,  region  identification 
table  and  initial  body  and  region  card 


CREW 

0 

MEMBERS 

22 

22 

5 3 SELL 
536 

15.20070 

9.00000 

-52.00000 

3.65000 

15.19957 

-52.00000 

7,15000 

AAA 

AAA 

537REC 

20.95614 

-52.00000 

-13.27500 

-11.50000 

.00000 

.00000 

AAC 

537 

-.00356 

-.00000 

11.07500 

-.00024 

4.00000 

-.00000 

AAC 

538RCC 

538 

15.20298 

2.00000 

-52.00000 

-3.45000 

-.00161 

-.00000 

5.00000 

AAB 

AAB 

539RCC 

539 

22.20721 

1.25000 

-52.00000 

-16.60000 

-.00463 

-.00000 

14.40000 

AAD 

AAD 

540RCC 

540 

8.20721 

1.25000 

-52.00000 

-16.60000 

-.00463 

-.00000 

14.40000 

AAE 

AAE 

541RCC 

541 

22.20673 

1.25000 

-50.75000 

-15.35000 

-.00073 

12.25000 

-.00000 

AAF 

AAF 

542RCC 

542 

8.20673 

1.25000 

-50.75000 

-15.35000 

-.00073 

12.25000 

-.00000 

AAG 

AAG 

543RCC 

543 

18.708889 

2.50000 

'-52.00000 

-21.85000 

-.00098 

16.50000 

-.00000 

AAH 

AAH 

544RCC 

544 

11.70889 

2.50000 

-52.00000 

-21.85000 

-.00098 

16.50000 

-.00000 

AAI 

AAI 

545RCC 

545 

18.70800 

1.50000 

-37.00000 

-21o85000 

.00482 

.00000 

-15.00000 

AAJ 

AAJ 

546RCC 

546 

11.70800 

1.50000 

-37.00000 

-21o85000 

.00482 

.00000 

-15.00000 

AAK 

AAK 

547ELL 

547 

0.00000 

9.00000 

-52.00000 

3.65000 

0.00000 

-52.00000 

7.15000 

AAA 

AAA 

548REC 

5o75000 

-52.00000 

-13.27500 

-11.50000 

.00000 

.00000 

AAC 

548 

0.00000 

-.00000 

11.07500 

0.00000 

4.00000 

-.00000 

AAC 

549RCC 

549 

0.00000 

2.00000 

-52.00000 

-3.45000 

0.00000 

-.00000 

5.00000 

AAB 

AAB 

550RCC 

550 

7.00000 

1.25000 

-52.00000 

-16.60000 

0.00000 

-.00000 

14.40000 

AAD 

AAD 

551RCC 

551 

-7.00000 

1.25000 

-52.00000 

-16.60000 

0.00000 

-.00000 

14.40000 

AAE 

Aj\E 

FIGURE  11:  Output  of  sample  problem  - title  card, 

card  and  body  table. 


parameter 


552RCC  7.00000  -50.75000 

552  1.25000 

553RCC  -7.00000  -50.75000 

553  1.25000 

554RCC  3.50000  -52.00000 

554  2.50000 

555RCC  -3.50000  -52.00000 

555  2.50000 

556RCC  3.50000  -37.00000 

556  1.50000 

557RCC  -3.50000  -37.00000 


557  1.50000 


536 

536 

0 

0 

537 

-537 

-536 

538 

538 

537 

0 

0 

539 

-537 

539 

0 

540 

-537 

540 

0 

541 

-539 

-537 

541 

542 

-540 

-537 

542 

543 

-537 

543 

0 

544 

-537 

544 

0 

545 

-543 

545 

0 

546 

-544 

546 

0 

547 

547 

0 

0 

548 

-548 

-547 

549 

549 

548 

0 

0 

550 

-548 

550 

0 

551 

-548 

551 

0 

552 

-550 

-548 

552 

553 

-551 

-548 

553 

554 

-548 

554 

0 

555 

-548 

555 

0 

556 

-554 

556 

0 

557 

-1 

-555 

557 

0 

FIGURE  11,  continued; 


15.35000 

OoOOOOO 

12.25000 

.00000 

AAF 

AAF 

15.35000 

0.00000 

12.25000 

.00000 

AAG 

AAG 

21.85000 

0.00000 

16.50000 

.00000 

AAH 

AAH 

21.85000 

0.00000 

16.50000 

* 

.00000 

AAI 

AAI 

21.85000 

0.00000 

.00000 

-15 

.00000 

AAJ 

AAJ 

21.85000 

0.00000 

.00000 

-15 

.00000 

AAK 

AAK 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

remainder  of  body  table  and  region  table 
from  output  of  sample  problem. 


•7) 


536 

10 

-0 

HEAD 

537 

11 

-0 

NECK 

538 

12 

-0 

BODY 

539 

13 

-0 

UPPER  ARM  RIGHT 

540 

14 

-0 

UPPER  ARM  LEFT 

541 

15 

-0 

LOWER  ARM  RIGHT 

542 

16 

-0 

LOWER  ARM  LEFT 

543 

17 

-0 

THIGH  RIGHT 

544 

18 

-0 

THIGH  LEFT 

545 

19 

-0 

LOWER  LEG  RIGHT 

546 

20 

-0 

LOWER  LEG  LEFT 

547 

10 

-0 

HEAD 

548 

11 

-0 

NECK 

549 

12 

-0 

BODY 

550 

13 

-0 

UPPER  ARM  RIGHT 

551 

14 

-0 

UPPER  ARM  LEFT 

552 

15 

-0 

LOWER  ARM  RIGHT 

553 

16 

-0 

LOWER  ARM  LEFT 

554 

17 

-0 

THIGH  RIGHT 

555 

18 

-0 

THIGH  LEFT 

556 

19 

-0 

LOWER  LEG  RIGHT 

557 

20 

-0 

LOWER  LEG  LEFT 

FIGURE  11,  continued:  region  identification  table  from  output 

of  sample  problem. 


NAMR(IOOO) 

JTYPE(2000) 

JBOD(30,30) 


NI  (30) 

LBOD(30,30) 

NM(30) 

2.  The  following  list  briefly  describes  the  sub- 
routines comprising  the  program  BRITL. 

BRITL  - is  the  main  program,  responsible  for  sequencing 
the  entire  calculation.  Most  of  the  input  is  read 
by  BRITL,  the  exceptions  being  the  thickness  data 
for  thick  bodies  (read  by  THICKN) , the  region  trans- 
formation data  (read  by  GENR) , and  the  pipe  data 
(read  by  PIPE) . As  the  input  body  and  region  tables 
are  processed,  BRITL  stores  the  output  bodies  and 
output  regions  on  TAPES  and  TAPES,  respectively » 

After  these  files  are  complete,  TAPES  and  TAPES  are 
copied  to  TAPE?  in  order.  The  input  region  identi- 
fication table  is  read  and  processed  and  the  results 
added  to  TAPE?.  BRITL  then  calls  the  subroutine 
STRIP  to  remove  unwanted  regions  (as  indicated  on  the 
region  identification  card) , to  remove  unused  bodies 
and  to  resequence  bodies  and  regions  consistently, 
starting  from  the  initial  body  number  and  initial 
region  number  given  on  the  last  input  card.  The 
final  output  is  on  TAPE?,, 


array  of  region  names  (3  characters)  versus 
region  nvimbers 

array  of  numbers  giving  the  body  type  and 
data  form  versus  body  number 

work  storage  area  for  Boolean  operations. 
OR'd  regions  are  stored  by  row.  The  columns 
of  a row  contain  bodies  which  intersect  to 
form  the  OR'd  region. 

number  of  intersecting  bodies  in  an  OR'd 
region  indexed  on  row  number  of  JBOD, 

similar  to  JBOD 

similar  to  NI  , 


PIPE (LI, L2)  - is  described  in  detail  in  III  D.  Reads  all 
but  the  leading  body  card  for  the  body  type  PIP,  pro- 
duces all  bodies  for  either  a thin  or  a thick  PIP, 
storing  the  data  in  locations  LI  to  L2  of  the  ASTER 
array. 

THICKN (Ll,L2 ,K)  - is  called  by  BRITL  when  an  input  body 
(whose  data  is  stored  in  MASTER  starting  at  K+1) 
is  determined  to  be  thick.  It  reads  the  thickness 
card  and  produces  the  new  body,  stores  the  data 
starting  at  LI  in  the  MASTER  array  and  returns  the 
final  location  L2  to  BRITL.  BRITL  then  calls 
WRITEB (Ll,L2)  to  print  and  punch  the  body  cards. 
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ARBT (MASTER (K+1) , IV, THICK)  - is  called  by  THICKN  to 

produce  the  second  body  for  a thick  ARB.  The  face 
descriptions,  contained  in  IV(4,6),  are  used  to 
identify  the  three  planes  whose  intersection  is  the 
vertex  under  consideration.  New  planes  are  obtained 
by  moving  each  plane  along  its  normal  without  ro- 
tation by  a distance  equal  to  the  thickness  for  the 
corresponding  facOo  The  new  vertex  is  the  inter- 
section of  the  three  new  planes o 

AQST  (THICK, LI, K,L2)  - is  called  by  THICKN  to  take  the 

input  AQS  data  (stored  in  MASTER  starting  at  K+1) , 
produce  a paraboloidal  surface  and,  using  the  local 
normal,  displace  each  of  the  original  points  along 
the  normal  by  the  input  thickness  to  produce  a second 
AQS  to  be  stored  in  MASTER  starting  at  LI  and  extend- 
ing to  L2,  the  latter  address  being  returned  to  THICKN. 


NOTR  - takes  a region  located  in  the  MASTER  array, 

negates  this  region,  leaves  the  results,  reduced 
to  the  conjunctive  form,  in  the  array  JBOD.  NOTR  sets 
the  arrays  (NI  (I)  ,I=1,N1),(  (JBOD  (J,  I)  ,J=1,NI  (I) ) ,I=1,N1) . 


RANDR  - intersects  two  regions,  each  in  the  conjunctive 
form,  to  produce  a final  region  in  the  conjunctive 
form.  One  of  the  original  regions  resides  in  the 
MASTER  array  and  the  other  is  stored  in  LBOD.  The 
results  are  ultimately  stored  in  LBOD,  NM  and  N2. 

DEROW  - is  called  by  NOTR  and  RANDR  to  complete  the 

simplification  of  a region  description.  In  essence, 
DEROW  removes  non-geometric  redundancies  by  removing 
OR'd  regions  which  are  explicitly  contained  in  other 
OR'd  regions. 


0RST0R(LB0D,N2,NM,MASTER(L1) ,L2)  - takes  the  contents  of 
the  array  ( (LBOD (L,M) ,L=1,NM (M) ,M=1,N2)  and  arranges 
it  in  MASTER  along  with  the  previously  stored  data 
for  the  current  region.  At  this  point  the  data  in 
LBOD  is  the  result  of  processing  a single  "OR"  of  an 
input  region  description.  THIS  "OR"  has  been  reduced 
(by  substitution  and  Boolean  manipulation)  to  N2  "OR"s, 
the  Mth  of  which  has  NM(M)  intersecting  bodies.  The 
NF  OR's  already  in  the  MASTER  array  are  stored  as  follows: 

MASTER (LI)  = NF 

MASTER(L1+M)  =NM(M),M=1, ,NF 


K-1 

MASTER (L1+N2+  Z NM(M)+I)= 
M=1 


body  number  (with  sign)  for  the  Ith 
body  of  the  Kth  "OR" . 
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ORSTOR  returns  a value  of  L2,  the  last  location 
occupied  in  the  MASTER  array. 

GENR(L1,L2)  - is  called  by  BRITL  to  produce  a new  region 
series  from  a model  region  series  by  transformation. 
GENR  reads  the  card  giving  the  number  of  regions  in 
the  series  and  the  cards  giving  the  4-point  trans- 
formation, deduces  a displacement  vector  and  a 
rotation  matrix  (or  a reflection  or  distortion  matrix) , 
and  then  applies  these  to  each  body  in  the  model 
regions  in  turn,  thus  producing  a new  sequence  of 
consecutively  numbered  bodies.  Concurrently,  the 
new  regions  are  constructed  by  substituting  new  body 
numbers  for  old  in  the  description  of  the  model 
regions.  GENR  calls  WRITES (LI , L2 ) for  each  new  body. 
The  region  write  is  invoked  by  GENR  with  a call  to 
WRITEF (NR, MASTER (LI) ,L2)  for  each  region  NR  in  the 
new  set. 

WRITES (Ll,L2)  - prints  and  punches  body  data  from  the 

storage  area  MASTER (Ll)  to  MASTER (L2).  The  punched 
(or  stored)  cards  are  in  the  format  acceptable  to  the 
GIFT  program. 

WRITEF (NR, MASTER (Ll) ,L2)  - takes  the  completed  region  des- 
cription for  region  NR  from  the  storage  area  MASTER (Ll) 
to  MASTER (L2)  and  constructs  the  region  cards  in  the 
format  acceptable  to  the  GIFT  program.  The  data  are 
both  printed  and  punched  (or  stored) . L2  is  calculated 
internally  and  returned  to  the  calling  routine. 


STRIP  - is  the  last  subprogram  called  by  BRITL  and  produces 
the  final  file  on  TAPE?.  STRIP  examines  the  region 
identification  cards  to  determine  which  regions  are  to 
be  retained.  Bodies  not  mentioned  in  these  regions  are 
deleted,  the  remaining  bodies  are  resequenced  and  the 
references  in  the  region  descriptions  are  corrected. 

The  region  cards  and  region  identification  cards  are  then 
resequenced.  Initial  body  number  and  initial  region 
number  are  obtained  from  the  last  input  card,  which  is 
read  by  STRIP. 

CROSS(A,B,C)  - computes  the  cross  product  of  B and  C and 
stores  the  answer  in  A. 

CUBIC (C,R,N)  - computes  the  roots  R of  the  cubic  equation 
whose  coefficients  are  input  as  C (unit  coefficient 
assumed  for  cubic  term) . Also  returns  the  number  of 
real  roots  N. 

DOT(A,B)  - computes  the  dot  product  of  A and  B. 
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SOLVE (A,B,N,KS>  - is  an  IBM  routine  that  solves  a set  of 
linear  equations,  order  N,  matrix  of  coefficients  A, 
constant  terms  in  vector  B,  The  solution  vector  is 
returned  in  B.  KS  is  a signal:  0 for  successful 
matrix  inversion,  1 for  failure, 

UNIT (A, B)  - normalizes  the  vector  A and  puts  the  results 
in  B without  destroying  A, 

UNPACK  - is  a utility  routine  used  in  connection  with  ARB 
data.  It  returns  an  index  related  to  the  number  of 
vertices. 
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III.  A COMGEOM  MODEL  OF  THE  M109A1  HOWITZER 

A COMGEOM  model  of  the  M109A1  Howitzer  has  been  prepared 
for  use  with  the  BRL-GIFT  program.  The  model  requires  939 
bodies  and  753  regions. 

Some  of  the  model  transformation  capabilities  of  BRITL 
(see  Section  II  of  this  report)  were  used  to  generate  copies 
of  repeated  structures  in  the  proper  position  and  orientation, 
e.g.  wheels  and  crew  members. 

A number  of  small  programs  were  written  for  the  checking 
and  correcting  of  body  parameters,  e.g.  mutual  orthogonality 
of  the  BOX  vectors  and  planar  faces  for  the  ARB. 

A major  checking  device,  namely  the  GIFT  program  itself, 
was  made  available  toward  the  end  of  the  contract  period,  and 
has  been  used  mainly  for  region  checking,  ioe.  to  discover 
overlaps. 
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DISTRIBUTION  LIST 


No.  or 

Organization  Copies  Organization 


No.  of 
Copies 

12  Connnander 

Defense  Dociamentation  Center 
ATTN;  DDC-TCA 
Cameron  Station 
Alexandria,  Virginia  22314 


1  Commander 

U.S.  Army  Materiel  Command 
ATTN;  AMCDL 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22333 

1 Commander 

U.S.  Army  Materiel  Command 
ATTN;  AMCRD,  BG  H.  A.  Griffith 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22333 

1 Commander 

U.S.  Army  Materiel  Command 
ATTN;  AMCRD-T 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22333 

1 Commander 

U.S.  Army  Aviation  Systems 
Command 

ATTN:  AMSAV-E 

12th  § Spruce  Streets 

St.  Louis,  Missouri  63166 

1 Director 

U.S.  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  California  94035 

1 Commander 

U.S.  Army  Electronics  Command 
ATTN:  AMSEL-RD 

Fort  Monmouth,  New  Jersey 
07703 


1 Director 

R&D  Tech  Support  Agency 
ATTN;  AMSEL-GG-E,  Mr.  C.  Goldy 
U.S.  Army  Electronics  Command 
Fort  Monmouth,  New  Jersey 
07703 

2 Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-R 

AMSMI-RDP,  Jerry  M. 

Brown 

Redstone  Arsenal,  Alabama 
35809 

2 Commander 

U.S.  Army  Tank  Automotive 
Command 

ATTN:  AMSTA-RHFL 

AMSTA-REMA, 

Mr.  Tasdemiroglu 
Warren,  Michigan  48090 

3 Commander 

U.S.  Army  Mobility  Equipment 
Research  § Development  Center 
ATTN-:  Tech  Docu  Cen,  Bldg.  315 

AMSME-RZT 

SMEFB-BS,  L.J.  Long 
Fort  Belvoir,  Virginia  22060 

1 Commander 

U.S.  Army  Armament  Command 
Rock  Island,  Illinois  61202 

1 Commander 

U.S.  Army  Picatinny  Arsenal 
ATTN:  SARPA-DW6,  Mr.  J.  Killen 

Dover,  New  Jersey  07801 

1 Commander 

U.S.  Army  Watervliet  Arsenal 
ATTN:  SARWV-RDD-SP 

J,  Wrzochalski 
Watervliet,  New  York  12189 
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DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


1 Commander 

U.S,  Naval  Missile  Center 
ATTN:  Code  5110,  B.E.  Nofrey 

Point  Mugu,  California  93042 

1 Commander 

U.S.  Naval  Weapons  Center 
ATTN:  Code  4565 

Mr.  C.  B.  Sandberg 
China  Lake,  California  93555 

1 Commander 

U.S.  Naval  Weapons  Laboratory 
ATTN:  Code  GAC,  Mr.  W.  Ralph 

Dahlgren,  Virginia  22448 

1 AFAL  (DLRV,  Mr.  J.  Flint) 
Eglin  AFB 
Florida  32542 

1 ASD  (XRHD,  G.B.  Bennett) 
Wright-Patterson  AFB 
Ohio  45433 


1 ASD  (ADDP,  Mr.  R.  Hilbrand) 
Wright-Patterson  AFB 
Ohio  45433 

1 ASD  (ENYW,  Mr.  J.  Howard) 
Wright-Patterson  AFB 
Ohio  45433 


Aberdeen  Proving  Ground 

Marine  Corps  Ln  Ofc 
Dir,  USAMSAA 

ATTN:  Dr.  J.  Sperrazza 

Mr.  L.  Bain 
Mr.  R.  Hoyt 
Mr.  R.  Marking 
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